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We extend the input-output formalism of quantum optics to analyze few-photon transport in waveguides 
with an embedded qubit. We provide explicit analytical derivations for one and two-photon scattering matrix 
elements based on operator equations in the Heisenberg picture. 
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I. INTRODUCTION 

In the context of quantum information technology, includ- 
ing quantum computing devices, understanding the interac- 
tion between a few-photon state and a two-level atom plays 
an important role [1—3]. The photons are a possible candi- 
date for the 'flying qubit' that carries the information, and the 
two-level atom constitutes the 'stationary qubit' where the fly- 
ing qubits are generated on demand and correlated with each 
other. 

Recently, there has been an increased activity in analyzing 
the properties of photons propagating in a waveguide coupled 
to a qubit — a two-level quantum mechanical system. Exper- 
imental demonstration of the control of single photons was 
made in a waveguide coupled to an optical cavity with an 
atom in its near field [4]. Similar effects were observed in the 
microwave domain, when low frequency photons in a trans- 
mission line were coupled to a superconducting qubit [5, 6], 
which later was shown to act as a photon amplifier [7]. 

To theoretically model such systems one needs to consider 
a continuous set of waveguide modes that are free to propa- 
gate in one dimension, either directly coupled to a multi-level 
system (referred to as an 'atom' in the paper), or indirectly 
coupled through an optical cavity with a discrete set of modes. 
Photon transport properties are non-trivial in these structures 
[8-11] which can be tailored to perform logic operations [12] 
or form a diode [13]. Exact solutions of one and two-photon 
scattering have first been reported in [9, 11]. 

The most widely used theoretical approach is to treat the 
set of equations in the Schrodinger picture, and apply the 
Lippmann-Schwinger formalism to calculate the reflection 
and transmission properties of the single and multi-photon 
states [10, 14-17]. An alternative technique is to use the re- 
duction formulas from field theory to calculate the scattering 
matrix of the system [18, 19]. Time-domain simulations that 
take the waveguide dispersion into account are also possible, 
and an interesting radiation trapping mechanism was recently 
predicted [20]. 

In this paper, we extend the input-output formalism [21, 22] 
of quantum optics — an Heisenberg picture approach origi- 
nally introduced to analyze the interaction between an atom 



in a cavity and a continuous set of electromagnetic states out- 
side of the atom-cavity system — to analyze the transport of 
few-photon states in a waveguide with an embedded qubit. In 
the input-output formalism one obtains a nonlinear set of op- 
erator equations based on the Hamiltonian of the system. For 
a coherent or a squeezed state input, this formalism has been 
extensively used to calculate various coherence properties of 
the output state of light. Here, we show that one can adopt this 
formalism to obtain exact results regarding one or two photon 
properties. To do so, we establish a relationship between the 
input-output formalism and the scattering matrix elements of 
the system. Our approach complements the existing theoreti- 
cal literature and bridges different analytical techniques. 

The paper is organized as follows: In Section II we intro- 
duce the Hamiltonian of the system. In Section III we build 
the link between the scattering theory and the input-output for- 
malism and continue in Section IV with the derivation of the 
one-photon transport properties. In Section V we show how 
to extend the calculations to the two-photon case. In Section 
VI we make observations on correlation function calculations 
based on coherent state inputs and end with our conclusions 
in Section VII. 



II. HAMILTONIAN 

We start by discussing the model Hamiltonian that we will 
use in this paper. As an illustration of the formalism, we con- 
sider a two-level atom coupled to a single polarization, single- 
mode waveguide [9], and treat the transport properties of few- 
photon states in such a system (Fig 1). The Hamiltonian, H, 
is defined as (ft = 1) 

H = H + H\ . 

Here Hq describes a chiral, i.e. one-way, waveguide where 
photons propagate in only one direction 
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and ctp, Zip are the annihilation and creation operators for the 
photons with a wavevector f3 respectively. In Appendix A we 
calculate the reflection and transmission probabilities for pho- 
tons in a waveguide where the fields propagate in both direc- 
tions and show that the results are straightforward extensions 
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FIG. 1. (Color online) Schematic representation of two photons in a 
waveguide, at frequencies k and p, moving to the right, towards a two 
level atom with energy levels \g) and \e). Q. is the separation between 
the energy levels. Coupling of the two level atom to the modes in the 
waveguide is proportional to t . Long horizontal lines denote the 
waveguide geometry. 



all these changes, we have 
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Throughout the paper, the labels for photon degrees of free- 
dom, for example k, p, refer to photon frequency. 



III. CONNECTION BETWEEN THE SCATTERING 
THEORY AND THE INPUT-OUTPUT FORMALISM 



of the chiral case. The operators obey the commutation rela- 
tion [hp, at,] = 6(J3 - /?')■ Hi describes the atom as well as the 
atom-waveguide interaction 

1 r 03 

Hi = + V J dy6 + ■ 

Here, Q is the atomic transition frequency, <x ± are the rais- 
ing and lowering operators for the two level atom and cr : = 
2cr + o-_ - I. V denotes the coupling strength between the 
atomic states and the waveguide modes. The derivation of 
the Hamiltonian is based on the dipole and the rotating wave 
approximations [23] as well as taking the continuum limit for 
field operators. The details of taking the continuum limit are 
discussed in Appendix B. 

It will be useful to have H in terms of the frequency of the 
photons instead of their wavevector, therefore, we linearize 
the waveguide dispersion around (J3o, ojq) as u>(J3) = coo+v g (J3- 
f3o) (see Fig 2). Notice that the total excitation operator 




commutes with H, i.e [ft, Ne] = 0. We could thus equiva- 
lent^ solve a system as described by 

H = H-a> Q N E =H Q + Hi (1) 

where 

d£v g 03-A,)a^ 

i r°° 

Hi = -Qo% + V d/3 (o- + a p + a^o-J) . 

Here £2 = £1 - o)o, and we also extended the lower limit of 
integration to -oo so that we can define the Fourier transform 
of operators in the next section. Since we will be dealing with 
states with wavevectors around fio, the extension of the inte- 
gration limit is well justified [24, 25]. Finally, we complete 
our transition from wavevectors to frequencies by defining 
<d = v g [5, and the operator a w = cip + p l JVg~, which satisfies 
the commutation relation [a w , a ,] = 6(a> - a>'). As a result of 



In a typical scattering experiment, various input states are 
prepared and sent towards a scattering region. After the scat- 
tering takes place, the outgoing states of the experiment are 
observed, and information about the interaction is deduced. 
The mathematical formulation of such a process is commonly 
made using the scattering matrix with elements of the form 

S PlP2 ,kik 2 = (PiP2\S\kik 2 ) 

where \k\k2) denotes the input states — here given as a two par- 
ticle state with energies (frequencies) ki and k 2 — and \pipz) 
the outgoing states. These input and output states are assumed 
to be free states in the interaction picture that exist long before, 
t — > -00, and long after, f — > 00, the interaction takes place. 
The S operator, then, is equal to the evolution operator in the 
interaction picture, Ui, from time -00 to +00, 

5= lim Ufa, t ) = lim e iHo ' l e" iH( ' l ^ o) e" iHo '° 

to— >-oo / ^-oo 
t\ — >oo /j — >oo 

where Hq is the non-interacting part of the Hamiltonian, and 
H = Hq +Hi is the total Hamiltonian. 1 In order to have a more 
compact notation, we will drop the limits and imply to — > -00 
and t\ — > 00. 

An equivalent way to describe the scattering is in terms of 
the scattering eigenstates \kik^) that evolve in the interaction 
picture from a free state either in the distant past or the distant 
future 

\k x k + 2 ) =U,{Q,to)\kik 2 ) = e iH '"e- iH "'"\kik 2 ) = Q + |Anfe> 
\k x k 2 ) sU^hmh) = e m e- iH ""\kik 2 ) = SlAhh). 

The interaction picture time evolution operators that relate 
scattering and free states are called the Miller wave opera- 
tors, Q. ± . The scattering operator can equivalently be written 
as S — Oif2 + . 2 It is also possible to write the scattering matrix 
elements as 

(pip 2 \S\kik 2 ) = (pip 2 \kik2). 



See [26] for more information about stationary scattering theory. [27] pro- 
vides a historical account of the developments related to the S -matrix. 
2 There is also an alternative definition of the scattering operator 5 ' = fi+nl 
which relates the incoming and outgoing scattering eigenstates, \k + ) = 
S'\k~), such that (p\S\k) = (p~\k + ) = (p _ |S'|r> = (p + \S'\k + ). See 
[28, 29] for details. 
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We should note that scattering eigenstates and the free states 
with the same quantum numbers have the same energies, that 
is H \kik 2 ) = E klk2 \k\k 2 ) and H\k\k±) = E k{kl \kikf) [26]. 

It is possible to denote the scattering matrix elements by an 
appropriate definition of input and output operators such that 



(pip 2 \kik%) = <0| flout(pi) a out (p 2 ) aUh) ahk 2 )\0) (4) 



where 




a m {k) = Cl + a k Q.' + 
a out (k) = Q.-a k £l'_ 



(6) 



have the property of creating input and output scattering eigen- 
states 

fl t(*)|0> = \k + ) 

«L(p)io> = i/o 

and the commutation relations 

[a in (k), a\ n (p)] = [a out {k),a' out (p)] = 8{k-p). 

We now relate the scattering theory, as briefly sketched 
above, to the input-output formalism [21, 22] of quantum op- 
tics. To do so, we start by recalling the definition of the input 
field operator [21] 



a m (t) 



i r 

V2^J ' 



' )e- ik( '- 

where a k (to) = e lHta a k e" m,<s is an operator in the Heisenberg 
picture. The relationship between a m (t) — which is defined 
in the input-output formalism — and a m (k) — which is defined 
above in (5) as a result of the scattering theory — can be deter- 
mined by noting that 



a m {t) 



= — f 

V2^ J 

= — f 

V2^ J 

=— r 

V2^ J 



dke iH '°a k e 



iHt -iH t „ JH a l -\Ht -ikt 



dke 



dk a in (k)e 



a k s 



(7) 



where in the second line we used the fact that [Hq, a k ] = -ka k 
to convert the a^e 1 *' term into e~ lH °'°a k e lHo ' . As a result, ai B (k) 
provides the spectral representation of a m {t) in the limit to — > 
-oo. Similarly, the output field operator in the input-output 
formalism 



flout(0 



V2^ J 



dka k (t\)e 



-m-h) 



is related to a out (k) in the scattering theory through 



a out(0 



- r 

V2^ J 



dk a out (k)e 



-ikt 



(8) 



in the limit fi — > oo. We have thus established a direct con- 
nection between the input-output formalism, and the scatter- 
ing theory. We should note that a different set of input and 



FIG. 2. (Color online) Linearization of a surface plasmon-like waveg- 
uide dispersion relation u)(fi) around a wavevector y8 is shown. The 
slope of the line is equal to the group velocity v g . The photon states 
in the text are assumed to have frequencies in the vicinity of a> so 
that the linearization is justified. 



output operators were defined in [30] with an aim to make a 
connection to correlation functions. In [31], a similar set of 
input-output operators were defined in order to relate two dif- 
ferent quantization schemes in dielectric media. To the best 
of our knowledge, the explicit link we provide above between 
the input-output formalism and the scattering theory has not 
been previously published in the literature. 



IV. SINGLE-PHOTON TRANSPORT 

Now that we know the relationship between the input- 
output formalism and the scattering theory, let us now calcu- 
late the S -matrix elements (p\S\k) between two single pho- 
ton states \k) and \p). Following the standard procedure, (see 
Appendix C), the input-output equations appropriate for the 
Hamiltonian in (1) are 



dJV 2 + 2 

— = -i \ ~(cr+ Gin - a! a-) - -N 
dt \ r ln r 



dcr_ 
~d7 



\l -cr z a u 



-—o~_- if2cr_ 

T 



CL'm —1 \\ —<T- 
T 



(9) 
(10) 
(ID 



where all operators are in the Heisenberg picture and hence 
they are all time-dependent. r~' = nV 2 /v g is proportional to 
the spontaneous emission rate. = cr + cr_ = (o~ z + l)/2 de- 
scribes the probability of having the atom in the excited state. 

The single-photon transport properties are described by the 
single photon S -matrix, which is related to the input and out- 
put operator by 



(p\S\k) = {0\a out (p)aUk)\0) 



~ f 



dt(Q\ a oat {t)\k + )e 



i/" 



where we used (8) to write a out (p) in terms of a O ut(0- It is 
therefore sufficient to first calculate (0| a O ut(0l^ + ) an d tnen 
perform an inverse Fourier transformation to determine the 
single-photon S -matrix. In the calculations to follow in this 
and the next section, we will go back and forth between 
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Fourier transforms of the operators, and we will explicitly use 
t, t' to imply time dependent operators and k\, 2 , p\, 2 to denote 
the time independent Fourier transformed pairs. 

The quantity (0| a out (t)\k + ) can be obtained by sandwiching 
(10) and (11) between the two states <0| and \k + ). We have 



£<0|o-_|* + > =iJ-(0\<r z a m \k+) - -<0|^_|^+> 
at V t t 



- iQ<0|cr_|/t + > 
\a oM \k + ) =<0| a in \k + ) 



T 



iicr_r>. 



Note that 



(0\a in (t)\k + ) 



I flm« <(*)!<>> 



-ikl 



V2^ 



by the use of (7) and 

(0\cr : a in (t)\k + ) 



\a in (t)\k + ) 



(12) 



(13) 



(14) 



(15) 



since |0) has an atomic part that is in the ground state. Using 
(14)— (15) in (12)— (13) results in a first order ordinary differen- 
tial equation. By solving it, we get 



<0|<r_|* + > = 
<0|fl out |fc + > = 



1 



-ikl 



V2^ (* - O) + i/r 

i ^ (fc-n)-i/T 

V2^ 6 (*-Q) + i/T 



(16) 
(17) 



After Fourier transforming (17), we obtain the single-photon 
S -matrix 



(p\S\k) = t k 6(k-p) 



(18) 



where 



tk = 



(fc - Q) - i/r 
(k - fij + i/r 



is the single-photon transmission coefficient. For subsequent 
calculations, we also define 



s k = 



2/T 



(Jfc - £2) + i/r 



that measures the excitation of the atom by the single-photon 
wave when normalized against the incident wave amplitude. 
tk and Sk are related by 



tk = 1 



Sk- 



These results for single-photon transport agree with [9, 11], 
where the scattering wavefunction was directly calculated 
through a real space formalism. 

The crucial step in the derivation above is (15) which takes 
advantage of the single-excitation nature of the input state. 
Formally, the same result can also be obtained by approxi- 
mately setting cr z = — 1 in (10), and thus linearizing the oper- 
ator equation. Such a procedure has been commonly adopted 



in many quantum optics calculations [32-34]. Typically, such 
an approximation is justified by assuming a so-called weak ex- 
citation limit, where the atom is assumed to be mostly in the 
ground state. Physically, in the case of single-photon trans- 
port, the weak excitation limit is valid, when a single-photon 
pulse has a duration that is much longer than the spontaneous 
lifetime of the atom. However, we emphasize that the weak- 
excitation limit is not always valid in general even for a single- 
photon pulse. It has been shown that for the Hamiltonian in 
(1), a single photon pulse with a duration comparable to the 
spontaneous emission lifetime can in fact completely invert an 
atom [35]. 

The formalism here removes the need for the assumption of 
weak-excitation limit when calculating single-photon proper- 
ties. In fact, we can directly calculate the excitation probabil- 
ity (k + \N\k + ) for the scattering eigenstate \k + ). N = cr+cr- and 
using (16) we have 

(k + \N\k + ) = (k + \o- + o--\k + ) = <*V + |0><0|o-_|Jfc + > 

= J_, |2 _ 1 2/r 

2n 2n(k-n) 2 + (l/T) 2 ' 

Here, we again have taken advantage of the fact that is a 
single-excitation state whereas <x + acting on any state except 
|0) would result in a multi-excitation state leading to a zero 
overlap with {k + \. More generally, we have 

(k + \o- z (t)\p + ) = {k%2o- + o-- - l)|p + > 

= 2<£V + |0><0k> + > - 5{k - p) 

= l e -«p-Qt s ls p -5(k-p) (19) 

71 

which will be useful when deriving the two-photon S -matrix. 

V. TWO-PHOTON TRANSPORT 

Our aim in this section is to calculate the two-photon 5- 
matrix based on the results we obtained for the single photon 
case. We first introduced the two photon S -matrix element in 
(4). We will begin by inserting an identity operator in between 
a ut(pi) and a oM (p 2 ) 

<0[ a out (pi) a out {p 2 ) a f m (ki) (& 2 )|0> 
= <0| floutOi) ( Jdk |£ + X*t + |) aUpi) al(h) al(k 2 )\0) 

and use the Fourier transform of (17) to simplify the expres- 
sion as 

= t pi (pi\ a out (p 2 ) fl^(fei) al n (k 2 )\Q). 
Using the Fourier transform of (1 1) we get 



Om(p2) -i-»/-0-_(p 2 ) 



aUh)aUk 2 )\G) 



= t Pi 6{p\ - ki)5(p 2 - k 2 ) + t Pi 5{p\ - k 2 )6(p 2 - h) 



-'^-t P ,{p\\o--(p 2 )\k l k + 1 ) 

T 
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where we used the orthogonality of the scattering eigenstates. 
Thus, to determine the two-photon S -matrix, we will need to 
calculate {py\o~-(i)\k\k%) and take its Fourier transform. 

Using (10), we obtain the differential equation that de- 
scribes (/?||cr_(?)|&i&J) 



-<ptl<r-(f)l*i**> 



(20) 



i v-0»ik«(«)flm(t)l*i^>-(- +in)<ptio-_(t)|ifciJ^>. 



If we can simplify the part that depends on cr. a in , we can then 
solve the differential equation. Since a\ n is an annihilation 
operator for scattering states, by using (7) we can write 

{p\\o- z {f)a m {t)\k x k1,} 
and then using (19) results in 



the algebraic Bethe ansatz [14], and the LSZ formalism in 
quantum field theory [18, 19]. The derivation here, however, 
is perhaps more elementary, and thus may serve to make such 
results more accessible. In addition, the results relate the pres- 
ence of the background fluorescence, to the excitation of the 
atoms. 



VI. COHERENT STATE COMPUTATION 

A traditional use of the input-output formalism is to calcu- 
late the correlation function when the input is in a coherent 
state. Here we briefly outline such a calculation for our sys- 
tem in order to contrast it with the single and two-photon cal- 
culations of the previous two sections. For this purpose, we 
consider a coherent input state |arjt), such that 

flin (fK> = ae"*K> 
and calculate, as an example, the G (1) correlation function 



(a A + |^ ut (Ofl out (Q|a A + ) 



V2^ 



V2^ 



which is what we were after. We can now solve the first order 
ordinary differential equation (20) in a way very similar to the 
derivation that led to (18). After some algebra and rearrange- 
ment we get 

- —^=~s kl+kl - Pl s* pi {s h + s k2 )e-' l(k ' +k ^' 
yln n 

+ -^=5{k 2 - Pl )s k] e-' lk " + - pi)s k ^ J . 

VZt V27T 



G m (t',t) 
Using (1 1), we have 

G il \t,t') =\a\ 2 &- m - f) + iae- [k < J-(cr + (t')) 



(22) 



where for any operator O, (O) = (at\0\at). 

Each of the expectation values in (22) can be calculated us- 
ing the input-output formalism. Taking the expectation values 
in (9) and (10) results in 



Taking the Fourier transform of the expression above gives us — (cr z (t)} — - i2 ^J—^ae lA '(cr + (f)) - a*e lk '(cr_(t))^ 



(ptlcr-ipimK) 

= --8{k\ +k 2 - p\ - p 2 )s P2 s* pj (s kl + s k2 ) 

+ s ki 5{k 2 - p\)5{ki - p 2 ) + s kl _6(ki - p\)5(k 2 - p 2 ). 
Lastly, using the relation t Pi s* p] = s Pl , we obtain 

<0[ a out (pi) a out (p 2 ) a\ a (ki) a\ n {k 2 )\0) 
= t kl t kl [6(k 2 - pi)5{ki - p 2 ) + 6(ki - pi)6(k 2 - p 2 )] 

1 [2 (2D 
+ l- a/ -6(ki +k 2 - pi- p 2 )s pi s P2 (s kl + %). 
n V r 

This final result agrees with previous calculations using ad- 
vanced techniques such as the Bethe ansatz 3 in real space [ 10], 



>c(0 + l> 



1 



' Equations (118)-(119) in [10] and equation (21) in this paper are the same 
with the following notational correspondence: T = 2/t, Ai = (k{ - fe)/2, 
A2 = (pi - P2)l% Ei = ki + k 2 , E 2 = pi + p 2 - 



-<cr + (f)) = I ill - - I <<r + (t)> - ia*e lk < ^ -(cr z (t)). 

Directly solving the equations above provides the values of 
<cr + (f')> and <cr_(r)) in (22), while the (cr+(f)cr-(t)) term can 
be computed using the quantum regression theorem. These 
calculations can be found in standard textbooks [22, 23], in 
sections related to the properties of resonance fluorescence, 
and we will not repeat them here. Instead, based on the out- 
line above, we make a few observations about the coherent 
state computations, as commonly done, and the one and two- 
photon computations as carried out in this paper. 

1 . The input-output formalism provides a set of nonlinear 
operator equations. Therefore, all computations, by neces- 
sity, involve the conversion of such operator equations into 
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ordinary differential equations for various operator matrix el- 
ements. While the coherent state computations typically in- 
volve taking expectation values in terms of the input states, 
the one and two-photon computations involve matrix elements 
that have different photon numbers. 

2. It is certainly reasonable to expect that the one or two- 
photon S matrices can be obtained by analyzing various cor- 
relation functions for a weak coherent state input. Indeed, the 
connection between the two-photon out wavefunction, and the 
g (2) correlation function, has been pointed out in [10] and it 
is likely that stronger connections exist. This will be carried 
out in future work. However, if the aim is to determine the 
S -matrix in the few -photon Fock state Hilbert space, the com- 
putation as discussed here should be far more direct. 

3. We emphasize that the few-photon computations yield 
the S -matrix in the few-photon Hilbert space, and thus pro- 
vide a complete description of all physical processes in the 
few-photon Fock state Hilbert space. In contrast, comput- 
ing G (1) or G (2) correlation functions alone do not completely 
specify the out state for a given incident coherent state in gen- 
eral. Certainly, in the majority of quantum optics experiments 
at present, one probes a quantum system with a coherent input 
state, and obtains information about the system by measuring 
different correlation functions. The coherent state computa- 
tions, as briefly reproduced above, are adequate to describe 
these experiments. However, these quantum systems are be- 
ginning to be considered as prospective devices which will 
eventually process quantum states [36, 37]. In such an engi- 
neering context, one ultimately has to be able to completely 
specify the output quantum states. It is in this respect that we 
hope the few-photon transport computations will prove to be 
valuable for future engineering applications. 



Appendix A: Two Mode Model 

In this section we will write the Hamiltonian for the case 
when photons are allowed to propagate in both directions 
within the waveguide. We will refer to this case as the two 
mode model. After introducing the Hamiltonian, we will use 
the results of Sections IV and V to calculate one and two 
photon reflection and transmission coefficients for right to left 
moving fields. 

When photons propagate to the right and to the feft, we will 
need to add extra terms to the Hamiltonian. We begin as we 
did in Section II and write 



Xoo j~0 



for the waveguide part of the Hamiltonian. The dispersion 
relation for the left moving modes a>c(J3) is the mirror image 
of the one for the right moving modes. We linearize the left 
and right branches of the dispersion relationship at f3 = +j6o to 
get oj r « o)q + v g (/3 - Pq) and cj t ~ (Oq - v g (/3 + fio). Follow- 
ing linearization, we extend the limits of integration to ±oo, 
make a change of variables /3 i-» /3 + /3q for the right and left 
waveguides respectively and define oj = v g /3, r w = rg +iSo / yvj, 

du u{rlr aJ -€ f Jo). (Al) 

oo 

The interaction part of the Hamiltonian is given by 

1 V f 00 

^ = -Qcr ; + — da> \cr + (r w + £J + (r+ + . 

(A2) 



Since the total excitation operator 



VII. CONCLUSION 

In this paper, we extend the input-output formalism of 
quantum optics to analyze one and two photon scattering in 
waveguides with a two-level atom inside. We develop the 
relationship between the input-output operators and the scat- 
tering theory which in turn enables us to analytically cal- 
culate the photon scattering matrix elements with minimum 
amount of algebra. We also contrast our calculations for few- 
photon Fock state transport with the conventional application 
of input-output formalism for coherent-state transport. This 
work helps us go beyond the correlation function analysis in 
input-output formalism, and leads to exact solutions for the 
scattering matrix elements. 



commutes with the Hamiltonian, we subtracted the term ojqNe 
from the Hamiltonian and set £1 = Q. - ojq in the derivation, 
mimicking the steps in Section II. 

Now that we have the Hamiltonian, we can write down the 
Heisenberg equations of motion and define the input-output 
operators for the fields as illustrated in detail for a chiral model 
in Appendix C. The equations for the annihilation operators 
are 



drjt) 
dt 

dt 



= -i[r w ,H] = -ia)r u - iVcr_ 
= -i[€ a ,H\ = +itoe a -iVo^ 
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where V — V/ ^vj. The definitions for the input and output 
operators for right going fields are the same as in Appendix C 
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r O ut(0 = r m (t) - i \\-o--(t) 
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where r is defined in (C6). Left going modes have a group 
velocity which is negative that of the right going modes and 
that leads to a negative sign in (Al). As a result, starting from 
the the definition of the input and output operators in (5)-(6), 
the input and output operators for left going modes have the 
form 



i-outCO = — f= fdwCOOe 
J2n J 



fcd(r-fi) 



V2^. 



where we note the change of sign in the frequency variable. 
Using these results we can show that 

dcr_ . [2 . [2 „ 2 

= 1 -»/ —cr, r m +1 -»/ — cr. t in — <x_ - ii2cr_ 

At V T V T T 

which is in a form similar to those that we get in temporal 
coupled mode theory [38]. 

We now have all the tools to solve for the scattering that 
takes place in the two mode model. Let us define even and 
odd combinations of the operators for the right and left propa- 
gating modes as 

*— a 



(even) a w = 



V2 ^ ' *" V2 
Using these definitions in (A1)-(A2) we can show 



(odd). (A3) 



H = j da> a'la [0 + °a^a 0J = H efi + H ofi 

~ = I ddJ (cr +fl(J + fl^cr_) = H e 
Vg J 



1 V2V 

Hi = -Qcr~ + 

2 



where we see that the interaction part of the Hamiltonian de- 
pends only on the even combination of modes. In Sections IV 
and V we solved for H = H s q + H e \ for a rescaled value of V. 
The odd part H D o is interaction free and hence is also solved. 
From (A3) we get 

flin/out(w) + flin/out(w) 



n n /out(w) = 



t (oS) 



V2 

fl ou t(-fr>) - aout(-a>) 
V2 



(A4) 



where we wrote the Fourier transforms of two-mode input- 
output operators in terms of the combinations of even and odd 
fields. 

The get the one photon reflection probability, we look at the 
scattering matrix element which corresponds to a right propa- 
gating input photon and a left propagating output photon 

<0K™t(p) a£(*)|o> 

= \(0\[a out (-p) ~ a out (-p)][al(k) + <£(*)]|0> 
= ^<0|aout(-p)<4(*)l°>- ^<0|flout(-p)4(*)l°> 
= ^(t k -l)S(p + k) = r k S(p + k). 



Here we used (A4) and (18) to get the one photon reflection 
coefficient f k . Similarly, the one photon transmission coeffi- 
cient t k is given by 

(0| r oM (p) rl(k)\0) = ife + \)5{p - k) = t k 6(p - k). 

Two photon calculations require adding another input- 
output pair. For instance, the scattering matrix element as- 
sociated with one photon scattering to the right, another to the 
left when two photons initially propagate to the right is given 
by 

(0|r out (/7i)4*(^2)^(A:i)^(fe2)|0> 

1 r + f 

= ^[(0| floutO'i) a out (-p 2 ) ajki) a in (k 2 )\0) 

- (0| flout(Pi) a ou t(-p 2 ) a^ik) aljk 2 )\0) 

- (0| flout(Pi) a out (-p 2 ) a^iki) al n (k 2 )\Q) 
+ (0| flout(Pi) a out (-p 2 ) al n (ki) al B (k 2 )\0) 
+ (0| flout(Pi) a out (-p2) a\ B {ki) a\ n {k 2 )\0) 

- (0| flout(pi) a out (-p 2 ) d^(fci) 0^(^2)10)] 

= i k fy 2 6(ki - p\)d(k 2 + p 2 ) + f k fy 2 5{k\ + p 2 )6(k 2 - p{) 

+ ^B8(k x +k 2 - pi + p 2 ) 
where from (21) 



1 2 

B = 1— -W —s Pl s- P2 (s k , + sic 2 ). 



We note that t' - t/2 due to an extra factor of V2 before V in 
the definition of Hi . These results agree with equations (52) 
and (130) in [10]. 



Appendix B: Hamiltonian in the continuum limit 

This section will summarize the steps taken to obtain the 
continuum form of the Hamiltonian from its discrete version. 
We will follow the approach in [24, 25]. 

The discrete variables are assumed to be for those in a one 
dimensional cavity of length L. The mode spacing in the cav- 
ity is given by AB = 2n/L. In this ID cavity, the free space 
electromagnetic Hamiltonian, Hq, is given by 



(Dp apdp 



with the commutator relationship [ap,a.p,] = 6pp. Now, 
we will convert the sum into an integral by the equivalence 
(A0 !/»)-»(/<#) to get 



H = ^ J dBioptfpC 
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The continuous mode operator tip is related to the discrete 
mode ap by 



Heisenberg equations of motion for the operators are 



hp = -J — cip which results in 
2n 



dySw^a' 



The commutator relationship [ap,al, ] = j^Spp in the limit 
L — > oo becomes 

To see this result, define /(0) = j^Spfi- Integrating /(0) will 
give 

/2n sr~t 2n L 

As a result, the correct Hamiltonian in the continuum limit is 
H = J d{3 oj(P) a^ap 

with [5g, at] = 5([3 - ft). It is then easy to show that 
where \/3) = a'^O), since 

In the discreet case 

#i = -Qxr z + — — Vfcr+Ofl + aj;Cr_) 

2 vi v " 

where V is the physical coupling constant. The factor Zr 1 / 2 
arises because the photon as created by has a normalization 
constant L~ l t 2 . In the continuum case we get 



1~ V L 2tt 

Hi =-fio-- + J — 

2 " tJl2tt y L 



J d/3(cr + ap + a^crJ) 



--Qcr z + , 
2 ~ V2lr 



y r 

V2^ J 



dp(cr + ap + a'„cr_). 



Thus, the coupling constants in the discrete (V) and the con- 
tinuum (V) cases differ by a factor of (2n) ■ 



Appendix C: Derivation of the input-output formalism 

Here we provide a derivation of the input-output equations 
(9)-(ll). This derivation closely follows [21, 22]. Based on 
the Hamiltonian (2)-(3), and the definition V = V/ Jv^, the 



dak 

1—7— = kak + Vct- 
dt 

dcr- 

i— ^- = Qcr_ - V I dkcr z ak 

i — - = 2V I dfc(-a!cr_ + cr+at). 
dt 



(CI) 
(C2) 

(C3) 



After multiplying (CI) by the integration factor exp(iAif), we 
integrate it from an initial time fo < t to get 



a k (t) = a k (t )e- lk ^ 
We define the input operator as 



-iV f dt'o--(t')e- ik{ '-'' ) . (C4) 

Jf 



flinW = f dfc fl ,(f )e- 1 * ( '-'° ) 



V2^ 

which satisfies the commutation relation 

[flin(0,<(^)] =<5(»-0- 

We further introduce a field operator 

O(0 = f dka k (t) 

V2tt J 

and integrate (C4) with respect to A: to get 

®(0 = flin(f) - i- V2mr_(t) = a in (t) - i J ^-(r-(f). ( C5 ) 
2 \ 2t 

Here, notice that we integrate over half the delta-function [21] 
which results in a factor of 1 /2 and r is defined as 

1 



= 7lV 2 . 



(C6) 



Furthermore, plugging (C5) into (C2) and (C3), results in 

dcr - • P 1 n 
= l -i/ —a, a m — cr_ — i£2cr_ 

df Vt t 

dN fl , 2 

— = -i y -(o-+ a m - a m o-_) - -at. 

Here iV = (<x, + l)/2. Thus the spontaneous emission rate 
is 2/t. We could have also directly calculated dN/dt from 
dcr_/df, since N = cr+cr-. 

Similarly, we integrate (CI) up to a final time t\ > t, and 
define an output operator 

fl O u,(0 = -L f djtflt^Oe-^) 
V2^ J 



which results in 



®(f) = a oat (t)+iJ—(r-(t). 

V 2t 

Combining (C5) and (C7), we finally obtain 



(C7) 



flout(f) = a iD (t) - i A / -cr„(f). 

T 
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